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Nuclear structure of 7 Be, 8 B and 7 ' 8 Li is studied within the ab initio no-core shell model (NCSM). 
Starting from high-precision nucleon-nucleon (NN) interactions, wave functions of 7 Be and 8 B bound 
states are obtained in basis spaces up to lOfifi and used to calculate channel cluster form factors 
(overlap integrals) of the 8 B ground state with 7 Be+p. Due to the use of the harmonic oscillator 
(HO) basis, the overlap integrals have incorrect asymptotic properties. We fix this problem in two 
alternative ways. First, by a Woods-Saxon (WS) potential solution fit to the interior of the NCSM 
overlap integrals. Second, by a direct matching with the Whittaker function. The corrected overlap 
integrals are then used for the 7 Be(p,7) 8 B S-factor calculation. We study the convergence of the 
S-factor with respect to the NCSM HO frequency and the model space size. Our S-factor results are 
in agreement with recent direct measurement data. We also test the spectroscopic factors and the 
corrected overlap integrals from the NCSM in describing the momentum distributions in knockout 
reactions with 8 B projectiles. A good agreement with the available experimental data is also found, 
attesting the overall consistency of the calculations. 

PACS numbers: 21.60.Cs, 21.30.Fe, 24.10.Cn, 25.40.Lw, 27.20.+n 



I. INTRODUCTION 

The Be(p,7) 8 B capture reaction serves as an impor- 
tant input for understanding the solar neutrino flux 0. 
Recent experiments determined the neutrino flux emit- 
ted from 8 B with a precision of 9% @. On the other 
hand, theoretical predictions have uncertainties of the 
order of 20% H Q. The theoretical neutrino flux de- 
pends on the 7 Be(p,7) 8 B S-factor that needs to be known 
with high precision. Many experimental and theoreti- 
cal investigations were devoted to this reaction. Exper- 
iments were performed using direct measurement tech- 
niques with proton beam and 7 Be targets |E IE as 
well as by indirect methods when a 8 B beam breaks up 
into 7 Be and proton in the Coulomb field of a heavy 
target [E • Theoretical calculations needed to extrap- 
olate the measured S-factor to the astrophysically rele- 
vant Gamow energy were performed with several meth- 
ods: the R-matrix parametrization [T(il |. the potential 
model Hdl Il2l Il3l | , and the microscopic cluster models 

0E1H 

In this work, we present the first calculation of the 
7 Be(p,7) 8 B S-factor starting from the ab initio wave func- 
tions of 8 B and 7 Be. We apply the ab initio no-core shell 
model (NCSM) 001. 

In this method, one considers 
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nucleons interacting by high-precision nuclcon-nuclcon 
(NN) potentials. There are no adjustable or fitted pa- 
rameters. Within the NCSM, we study the binding en- 
ergies and other nuclear structure properties of Be, 8 B 
as well as 7 ' 8 Li, and calculate overlap integrals for the 8 B 
and 7 Be bound states. Due to the use of the harmonic- 
oscillator (HO) basis, we have to correct the asymptotic 
behavior of the NCSM overlap integrals. This is done 
in two alternative ways. First, by fitting Woods-Saxon 
(WS) potential solutions to the interior part of the NCSM 
overlap integrals under the constraint that the experi- 
mental 7 Be+p threshold is reproduced. Second, by a 
direct matching of the NCSM overlap integrals and the 
Whittaker function. The corrected overlap integrals are 
then utilized to calculate the 7 Be(p,7) 8 B S-factor as well 
as momentum distributions in stripping reactions. We 
pay special attention to the convergence of the S-factor 
with respect to the NCSM model space size and the HO 
frequency. 

In Sect. UH we present the NCSM results for 7 Be, 8 B 
and 7,8 Li energies, ground-state radii and electromag- 
netic moments and transitions. The calculation of clus- 
ter form factors and the correction of their asymptotics 
is described in Sect. II I II A test of the corrected NCSM 
overlaps and spectroscopic factors for the momentum dis- 
tributions in the stripping reaction ( 7 Bc + p) + A — > 
7 Be + X is discussed in Sect. IIVI The S-factor calcula- 
tion with sensitivity and convergence studies is presented 
in Sect. Conclusions are drawn in Sect. IVll 
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II. NUCLEAR STRUCTURE OF 7 BE, 7 LI, 8 B 
AND 8 LI 

A. Ab initio no-core shell model 

In the NCSM, we consider a system of A point-like 
non-relativistic nucleons that interact by realistic two- 
or two- plus three-nucleon interactions. The calculations 
are performed using a finite harmonic oscillator (HO) ba- 
sis. As in the present application we aim at describing 
loosely bound states, it is desirable to include as many 
terms as possible in the expansion of the total wave func- 
tion. By restricting our study to two-nucleon (NN) inter- 
actions, even though the NCSM allows for the inclusion 
of three-body forces [l^, we are able to maximize the 
model space and to better observe the convergence of our 
results. The NCSM theory was outlined in many papers. 
Here, we only briefly review the main points for the case 
when the Hamiltonian is restricted to just a two-nucleon 
interaction. 

We start from the intrinsic two-body Hamiltonian for 
the A-nucleon system Ha = T re i+V, where T re i is the rel- 
ative kinetic energy and V is the sum of two-body nuclear 
and Coulomb interactions. Since we solve the many-body 
problem in a finite HO basis space, it is necessary that we 
derive a model-space dependent effective Hamiltonian. 
For this purpose, we perform a unitary transformation 
(17L \18L l20l l21j of the Hamiltonian, which accommodates 
the short-range correlations. In general, the transformed 
Hamiltonian is an A-body operator. Our simplest, yet 
non-trivial, approximation is to develop a two-particle 
cluster effective Hamiltonian, while the next improve- 
ment is to include three-particle clusters, and so on. The 
effective interaction is then obtained from the decoupling 
condition between the model space and the excluded 
space for the two-nucleon transformed Hamiltonian. De- 
tails of the procedure are described in Refs. dill. The 
resulting two-body effective Hamiltonian depends on the 
nucleon number A, the HO frequency f2, and N^ax, the 
maximum many-body HO excitation energy defining the 
model space. It follows that the effective interaction ap- 
proaches the starting bare interaction for iV max — > 00. 
Our effective interaction is translationally invariant. A 
significant consequence of this fact is the factorization 
of our wave functions into a product of a center-of-mass 
|M7 component times an internal component. 

Our most significant approximation here is the use 
of the two-body cluster approximation to the effective 
many-body Hamiltonian. Our method is not variational 
so higher-order terms may contribute with either sign to 
total binding. Hence, evaluating the dependence on the 
basis-space parameters help calibrate our convergence. 

Once the effective interaction is derived, we diagonal- 
ize the effective Hamiltonian in a Slater determinant HO 
basis that spans a complete N max hH, space. This is a 
highly non-trivial problem due to very large dimensions 
we encounter. To solve this problem we have used a spe- 
cialized version of the shell model code Antoine 1221 l23l , 



recently adapted to the NCSM |2l|2^. This code works 
in the M scheme for basis states, and uses the Lanczos al- 
gorithm for diagonalization. Its basic idea is to write the 
basis states as a product of two Slater determinants, a 
proton one and a neutron one. Matrix elements of opera- 
tors are calculated for each separate subspace (one-body 
for the proton-neutron, two-body for the proton-proton 
and neutron- neutron). The performance of the code is 
the best when the ratio between the number of proton 
plus neutron Slater determinants and the dimension of 
the matrix is the least. It happens when the number of 
proton Slater determinants is equal the number of the 
neutron Slater determinants. The number of iterations 
needed to converge the first eigenstates is significantly 
reduced by the implementation of a sophisticated strat- 
egy for selecting the pivot vectors. The difficulty of the 
no-core calculations (in which all nucleons are active) is 
that the number of shells and, consequently, the number 
of matrix elements that are prccalculated, becomes very 
large. One has to handle a huge number of operators. 
This is the reason why it has been necessary to write a 
specialized version of the code. 

A recent development of the NCSM is the ability to 
further process the wave functions, resulting from the 
shell-model calculation, to obtain channel cluster form 
factors [2(j and, consequently, the spectroscopic factors. 



B. Nucleon-nucleon interactions 

Two different, high-precision (i.e. such that provide 
a perfect fit to two-nucloen data) NN interactions have 
been used in this study: the CD-Bonn 200 |27| and the 
INOY (Inside Nonlocal Outside Yukawa) |2S U 29 j poten- 
tials. 

The CD-Bonn 2000 potential as well as its earlier 
version [3(j is a charge-dependent NN interaction based 
on one-boson exchange. It is described in terms of co- 
variant Fcynman amplitudes, which are non-local. Con- 
sequently, the off-shell behavior of the CD-Bonn inter- 
action differs from local potentials which leads to larger 
binding energies in nuclear few-body systems. 

A new type of interaction, which respects the local 
behavior of traditional NN interactions at longer ranges 
but exhibits a non-locality at shorter distances, was re- 
cently proposed by Dolcschall et al. [2^, [2i|. The au- 
thors explore the extent to which effects of multi-nucleon 
forces can be absorbed by non-local terms in the NN 
interaction. Their goal was to investigate if it is pos- 
sible to introduce non-locality in the NN interaction so 
that it correctly describes the three-nucleon bound states 
3 H and 3 He, while still reproducing NN scattering data 
with high precision. The so called IS version of this in- 
teraction, introduced in Rcf. [2^], contains short-range 
non-local potentials in 1 S'o and 3 Si — 3 Di partial waves 
while higher partial waves are taken from Argonne Dig . In 
this study we are using the IS-M version, which includes 
non-local potentials also in the P and D waves |2|j. It 
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is important to note that, for this particular version, the 
on-shell properties of the triplet P-wave interactions have 
been modified in order to improve the description of 3N 
analyzing powers. The 3 Pq interaction was adjusted to 
become less attractive, the 3 Pi became more repulsive, 
and the 3 P2 more attractive. Unfortunately, this gives a 
slightly worse fit to the Nijmegen 3 P phase shifts. 



C. 7 Be and 7 Li 

Our calculations for both A = 7 and A = 8 nuclei were 
performed in model spaces up to 10M7 for a wide range 
of HO frequencies. We then selected the optimal HO fre- 
quency corresponding to the ground-state energy mini- 
mum in the 10M7 space and performed a 12M7 calculation 
to obtain the ground-state energy and the point-nucleon 
root-mean-square (rms) radii. The overlap integrals as 
well as other observables were, however, calculated only 
using wave functions from up to 10M7 spaces. 

The 7 Be ground-state energy dependence on the HO 
frequency for different model spaces is shown in Figs. ^ 
and for the CD-Bonn 2000 and the INOY NN poten- 
tials, respectively. We observe a quite different conver- 
gence trend for the two potentials. For the INOY, the 
convergence is very uniform with respect to the HO fre- 
quency with substantial changes with N max . The con- 
vergence with increasing N max is evident as also seen in 
the inset of the Fig. [5] It is straightforward to extrap- 
olate that the converged INOY ground-state energy will 
slightly ovcrbind 7 Bc. The ground-state energy conver- 
gence for the CD-Bonn is quite different with a stronger 
dependence on the frequency, with minima shifting to 
smaller frequency with basis size increase, and an overall 
weaker dependence on iV max . Contrary to the INOY, the 
CD-Bonn undcrbinds 7 Bc by more than 3 MeV, which is 
typical for the standard high- precision NN potentials. 

From these results wc select the optimal frequency 
HQ = 12 MeV for the CD-Bonn 2000 and hfl = 16 MeV 
for the INOY potential. Corresponding spectra obtained 
using these frequencies are then shown in Figs. [3] and ^ 
The energies, radii and electromagnetic observables are 
summarized in Tables [I] and [HI where we also include 
the 7 Li results. We obtain the same level ordering for 
7 Be and 7 Li which is also the same for both NN poten- 
tials with the exception of a reversal of the 7/2^ and 
1/2^ levels. Our CD-Bonn ordering is in agreement with 
experiment for the 9 lowest levels in 7 Li. In 7 Be, the ex- 
perimental 7/2^ and 3/2^~ levels are reversed compared 
to our results and the situation in 7 Li. Our calculated 
spectra show a good convergence with the basis size en- 
largement as well as a good stability with respect to the 
frequency change. An interesting feature is a different 
level spacing of the 7/2j~, 5/2^ and 5/2^ levels for the 
CD-Bonn and INOY calculations. The INOY results are 
similar to those obtained when a three-body interaction 
is included in the Hamiltonian [l9l | . This is not surprising 
as discussed in Subsect. Ill Bl 



Concerning the magnetic moments and Ml transitions, 
we can see a very little dependence of the calculated val- 
ues on the HO frequency or the basis size. Also the 
two NN potentials give very similar results that are in a 
reasonable agreement with experiment. Concerning the 
radii and quadrupole moments, the calculated values in 
general increase with increasing basis size and decreasing 
frequency. This is in part a consequence of the incorrect 
asymptotics of the HO basis. The fastest convergence 
for the radii and quadrupole moment occurs at a smaller 
HO frequency. In the CD-Bonn calculations for Ml = 11 
and 12 MeV, the radii are close to experimental values. 
The quadrupole moments are still underestimated. For 
the INOY potential we observe underestimation of both 
radii and quadrupole moments. This is not unexpected 
as the INOY NN potential also underpredicts the 4 He 
radius. 



D. 8 B and 8 Li 

Our 8 B ground-state energy results are presented in 
Figs. [5]and[f)]as well as in Table ITTT1 where we also show 
the excitation energies and the corresponding 8 Li results. 
The basis size dependence of the 8 B spectra calculated 
using the optimal HO frequencies is shown in Figs. [7| 
and GO for the CD-Bonn and INOY NN potentials, re- 
spectively. Similar conclusions can be drawn as for the 
A = 7 nuclei concerning convergence. The INOY NN 
potential gives binding energy close to the experimen- 
tal one and a small overbinding is expected based on 
extrapolation of our N max dependence. The CD-Bonn 
undcrbinds both 8 B and 8 Li by about 5 MeV. We note 
that the A = 8 nuclei, with emphasis on 8 Be, were also 
investigated within the NCSM in Ref. 0. The present 
8 B and 8 Li CD-Bonn binding energy results are basically 
identical to those of Ref. [24| with small differences due 
to the use of different versions, j2]j vs. [3(| , of the CD- 
Bonn NN potential. It is interesting to note that both 
employed NN potentials predict 8 B unbound, contrary to 
experiment. The problem is less severe for the INOY NN 
potential. This suggests that a three-nucleon interaction 
is essential to reproduce the experimental threshold. This 
may appear as a problem in view of the present applica- 
tion to 7 Be(p,7) 8 B reaction. However, since our basis 
has incorrect asymptotics in the first place, we make use 
of only the interior part of our ab initio wave functions 
that are presumably unaffected by the wrong threshold 
placements. This is discussed in the next section. 

Concerning the ground-state energy systematics of the 
A = 7 and A = 8 nuclei, wc note that recent Green's 
function Monte Carlo (GFMC) calculations for 7 Li and 
8 Li using the Argonne potentials |3^| found qualitatively 
similar differences between the calculated and experimen- 
tal values as we observe here. 

Concerning the excitation energies, a noticeable differ- 
ence between the CD-Bonn and the INOY predictions for 
the low-lying levels is a reversed order of the 0^ and 3~/~ 



4 



states. We note, however, that the state has not been 
observed experimentally. 

As seen from Table Hvl the radii and quadrupole mo- 
ments are substantially larger and closer to experiment 
in our CD-Bonn calculations. Contrary to the 7 Be- 7 Li 
case, here we observe an interesting difference between 
the two NN potentials for the magnetic moment predic- 
tion. The INOY NN potential gives the magnetic mo- 
ment of 8 Li greater than that of 8 B in agreement with 
experiment, while the CD-Bonn predicts the opposite. 
Interestingly, almost identical 8 Li and 8 B magnetic mo- 
ments as we obtained using the CD-Bonn were reported 
in Ref. ^(| calculated within a cluster model applied to 
the 7 Be(p,7) 8 B reaction. Clearly, our results suggest that 
the A = 8 magnetic moments are sensitive to a presence 
of thrcc-nuclcon interaction in the Hamiltonian (that is 
mocked up to certain degree by the INOY NN potential) . 
As to the B(M1;1 + — > 2 + ) transition, both potentials 
gives about the same result that is in agreement with the 
8 Li experimental value but by almost a factor of three 
smaller than the experimental value for 8 B. We note the 
large experimental error bar of the latter. It appears that 
the 8 B experimental values is inconsistent with both the 
analogous 8 Li value and our calculations. 



III. CLUSTER FORM FACTORS FROM AB 
INITIO WAVE FUNCTIONS 



Overlap functions and spectroscopic factors 
obtained in the model space 



Detailed knowledge of nuclear structure is important 
for the description of low-energy nuclear reactions. As 
the first step in the application of the NCSM to low- 
energy nuclear reactions, one needs to understand the 
cluster structure of the eigenstates. This is achieved by 
calculating the channel cluster form factors (or overlap 
integrals, overlap functions). The formalism for calcu- 
lating the channel cluster form factors from the NCSM 
wave functions was developed in Ref. |26| . Here we just 
briefly repeat and adapt a part of the formalism relevant 
for our present application. 

We consider a composite system of A nucleons, i.e. 8 B, 
a miclcon projectile, here a proton, and an A — 1-nuclcon 
target, i.e. 7 Be. Both nuclei are assumed to be described 
by eigenstates of the NCSM effective Hamiltonians ex- 
panded in the HO basis with identical HO frequency and 
the same (for the eigenstates of the same parity) or dif- 
fering by one unit of the HO excitation (for the eigen- 
states of opposite parity) definitions of the model space. 
The target and the composite system is assumed to be 
described by wave functions expanded in Slater determi- 
nant single-particle HO basis (that is obtained from a 
calculation using a shell model code like Antoine). 



Let us introduce a projectile-target wave function 
(6... U-*r A®":^ ;S r ) 



; (jmh Mr | JM) (Zm ; \ m 8 \jm) s -^l 
x Y lmi (f) Xms (H • ■ ■ U-M - lahAh) , (1) 



where (£i . . . £a-2|^4- — \aI\M\) and Xm 3 are the tar- 
get and the nucleon wave function, respectively. Here, 
/ is the channel relative orbital angular momen- 
tum, £ are the target Jacobi coordinates and r = 
-fi^j (fi + r*2 + . . . + ta-i) — Ta describes the relative 

distance between the nucleon and the center of mass of 
the target. The spin and isospin coordinates were omit- 
ted for simplicity. 

The channel cluster form factor is then defined by 

^-w 1 « = (^AJ|^j^;<5 r > ) (2) 

with A the antisymmetrizer and |^4AJ) an eigenstate of 
the A-nucleon composite system (here 8 B). It can be 
calculated from the NCSM eigenstates obtained in the 
Slater-determinant basis from a reduced matrix clement 
of the creation operator: 



(nlOQl\QOnll) i J 

X SDMAJ||4yp-laJi>SD • (3) 

The subscript SD refers to the fact that these states 
were obtained in the Slater determinant basis, i.e. by 
using a shell model code, and, consequently, contain spu- 
rious center-of-mass components. A general HO bracket, 
which value is simply given by 



(nl00l\00nll)_L- = (-1)' 



A-l 
A 



(4) 



then appears in Eq. J3J in order to remove these compo- 
nents. The R n i (r) in Eq. © is the radial HO wave func- 
tion with the oscillator length parameter b = A-i m ^ , 

where m is the nucleon mass. 

A conventional spectroscopic factor is obtained by in- 
tegrating the square of the cluster form factor: 



A\J 



(5) 



We calculated the ( 8 B| 7 Be+p) channel cluster form fac- 
tors and the spectroscopic factors for the bound states 
of 8 B and 7 Be from the NCSM wave functions obtained 
in model spaces up to N mgx = 10 (10M7) in a wide HO 
frequency range and using both the CD-Bonn and the 
INOY NN potentials. The most important channels are 
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the p- waves, 1 = 1, with the proton in the j = 3/2 and 
j = 1/2 states. Our selected spectroscopic factor results 
are summarized in Table [3 It should be noted a very 
weak dependence of the spectroscopic factors on either 
the basis size or the HO frequency or the NN potential. 

We note that the ( 8 Li| 7 Li+n) spectroscopic fac- 
tors were obtained from 7 Li(d,p) 8 Li stripping measure- 
ments through a distorted-wave Born approximation 
(DWBA) analysis 0IU. The NCSM ( 8 Li| 7 Li+n) and 
( 8 B| 7 Bc+p) spectroscopic factors are almost identical 
and overestimate those extracted in Ref. |34l |3f| with 
the ground-state values of about 0.9. At the same time, 
the NCSM spectroscopic factors are in good agreement 
with those obtained in R-matrix analysis by Barker |lOj| . 
Also, they are consistent with the spectroscopic factors 
from the microscopic three-cluster model lift as well as 
with the Cohen-Kurath spectroscopic factors |3(j with no 
center-of-mass correction. 



B. Correction of the overlap- function asymptotics 

The dominant p-wave j = 3/2 overlap integral for the 
8 B and 7 Be ground states obtained using the CD-Bonn 
NN potential in the 10M7 model space and the HO fre- 
quency of Ml = 12 MeV is presented in Fig. [5] by the 
full line. Despite the fact that a very large basis was 
employed in the present calculation, it is apparent that 
the overlap integral is nearly zero at about 10 fm. This 
is a consequence of the HO basis asymptotics. A sepa- 
rate issue is the incorrect threshold placements obtained 
from the high precision NN potentials that we employ. 
The proton capture on 7 Be to the weakly bound ground 
state of 8 B associated dominantly by the El radiation is 
a peripheral process. Consequently, the overlap integral 
with an incorrect asymptotic behavior cannot be used to 
calculate the S-f actor. 

It is our expectation, however, that the interior part 
of the overlap integral as obtained from our large-basis 
NCSM calculation is realistic. This is supported by the 
fact that the overlap integrals obtained using the CD- 
Bonn 2000 and the INOY NN potentials are similar de- 
spite the fact that the two potentials predict different 
binding energies and thresholds. Further, we calculated 
overlap integrals using a Hamiltonian that included the 
three-nucleon interaction. In those calculations, for tech- 
nical reasons limited to small model spaces up to 6?iQ, 
we obtained binding energies closer to experiment. Still, 
the interior part of the overlap integrals was almost the 
same as that calculated here from the CD-Bonn 2000 in 
the same model space and identical HO frequency. It is 
then straightforward to correct the asymptotic behavior 
of the overlap integral and fix in this way both the HO 
basis and the wrong threshold placement issues at the 
same time. One possibility wc explored utilizes solutions 
of a Woods-Saxon potential. In particular, we performed 
a least-square fit of a WS potential solution to the inte- 
rior of the NCSM overlap in the range of — 4 fm. The 



WS potential parameters were varied in the fit under 
the constraint that the experimental separation energy 
of 7 Be+p, i.e. 0.137 MeV, was reproduced. In this way 
we obtain a perfect fit to the interior of the overlap inte- 
gral and a correct asymptotic behavior at the same time. 
The resulting overlap function is presented in Fig. HJby 
the dashed line. 

Another possibility is a direct matching of logarithmic 
derivatives of the NCSM overlap integral and the Whit- 
taker function: 

^-ln(r 9lj (r)) = ^ln(C l3 W^ l+1/2 (2k r)) , (6) 

where rj is the Sommerfeld parameter, fco = \/2^EQ/h 
with fi the reduced mass and Eq the separation energy 
(here E = 0.137 MeV). The NCSM overlap integral 
gij(r) is defined in Eq. |(2J), and, according to Eq. J2J), 
it is expanded in terms of the HO radial wave functions, 
i.e. rgi 3 (r) = J2 n s n ijTl n i(r) withK n i{r) = rR n i(r). The 
spectroscopic factor (JEJl is then given by Sij = J2 n s nij ■ 
For simplicity, we suppressed all quantum numbers ex- 
cept the orbital momentum I, the total proton angular 
momentum j and the radial HO quantum number n. 
With the help of 

d ^ , s f 2n + l + l r \ . , 
ar \ r o z / 

-^n(n + l + i)Vi,(r) , (7) 

it is straightforward to calculate the derivative of r gij (r) . 
Since asymptotic normalization constant (ANC) Cij can- 
cels out, there is a unique solution of Eq. © at r = R m . 
For the discussed overlap presented in Fig. |5J we found 
R m = 4.05 fm. Finally, by matching the value of the 
NCSM overlap and the Whittakcr function at R m , we de- 
termine the ANC. The corrected overlap using the Whit- 
taker function matching is shown in Fig. [5] by a dotted 
line. In general, we observe that the approach using the 
WS fit leads to deviations from the original NCSM over- 
lap starting at a smaller radius. In addition, the WS 
solution fit introduces an intermediate range from about 
4 fm to about 6 fm, where the corrected overlap deviates 
from both the original NCSM overlap and the Whittakcr 
function. Perhaps, this is a more realistic situation com- 
pared to the direct Whittaker function matching. In any 
case, by considering the two alternative procedures we 
are in a better position to estimate uncertainties in our 
S-factor results. This is discussed later. 

In the end, we re-scale the corrected overlap functions 
to preserve the original NCSM spectroscopic factor (given 
in Table lv|) . In general, we observe a faster convergence 
of the spectroscopic factors than that of the overlap func- 
tions. With increasing basis size (iV max ), the maximum 
of the overlap functions decreases while the tail extends 
with the integral of the square approximately conserved. 
This is demonstrated in Fig. 1101 see further discussion 
below. The corrected overlap function should represent 
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the infinite space result. By re-scaling a corrected over- 
lap function obtained at a finite N max , we approach faster 
the infinite space result. 

The same procedure is applied to other relevant chan- 
nels. In Fig. ^2 we present the NCSM overlap inte- 
gral and its corrected form for the other p-wave channel 
with proton in j = , and 8 B and 7 Be in their ground 
states. Obviously, the WS parameters as well as the 
matching radii R m obtained for the two channels are 
different (in the case of the present j = h channel we 



found R m = 4.23 fm). In Table rVTI we show the fitted 
WS potential parameters obtained in the two discussed 
cases together with parameters corresponding to a p-wave 
channel with 7 Be in the first excited state. In this last 
case, the separation energy is E = 0.57 McV. We use the 
definition of the WS potential as given, e.g. in Eqs. (5-7) 
of Ref. j3]J . Typically, the central potential parameters 
Ro, ao are well constrained in the fit, while the spin-orbit 
potential parameters are obtained with some uncertainty. 
Their values then exhibit more variation from channel to 
channel. The strength of the central potential Vq was 
re-adjusted at every step during the fit to reproduce the 
experimental separation energy. We note that parameter 
values in Table IVII are rounded compared to the results 
from our fitting procedure. One needs to fine-tune, e.g. 
the Vb, to reproduce accurately the respective separation 
energy in order to use the parameters independently. 

It should be noted that the Woods-Saxon potential 
here is just a tool to represent the interior part of the 
NCSM overlap integral as accurately as possible and cor- 
rect its asymptotic form at the same time. The range 
used in the least-square fit is not arbitrary and varies 
from channel to channel. The aim is to use as large range 
as possible, while at the same time preserve the NCSM 
overlap integral as accurately as possible in that range. 
A combination of eye-guided evaluation with a quantita- 
tive condition of minimizing the least-square function per 
fitted point helps to determine the largest possible range. 
Finally, let us repeat that the alternative procedure of the 
direct Whittakcr function matching is completely unique. 

The corrected overlap integrals then serve as the input 
for the momentum distribution and the S-factor calcula- 
tions as described in the following sections. 

The basis size dependence of the overlap integrals both 
the original NCSM and the corrected ones (using the WS 
solution fit procedure) can be judged from Fig. 1101 where 
we compare results obtained in model spaces from &MI to 
10MI using the CD-Bonn 2000 NN potential and a fixed 
HO frequency of Ml = 12 MeV. With an increase of the 
basis, the maximum of the overlap decreases and its tail 
extends to larger distances. The corresponding spectro- 
scopic factors differ very little, see Table Ivl The change 
from &MI to \QMl is smaller than that between QMl and 
8M1, a sign of convergence. In Fig. ^1 wc then com- 
pare the overlap integrals, both original and corrected, 
obtained in the \QMl model space using three different 
HO frequencies. With a decrease of the HO frequency, 
the maximum of the overlap decreases and shifts to a 



larger distance. At the same time, the tail extends to 
larger distances as well. Finally, in Fig. ^] we compare 
overlap integrals obtained using the CD-Bonn 2000 and 
the INOY NN potentials in the same lOMl model spaces 
and the same HO frequency of Ml — 14 MeV. While 
the spectroscopic factors are almost identical, the CD- 
Bonn overlap integral has a smaller maximum shifted 
to a larger distance and also extends further. The pa- 
rameters of the WS potential that produce the INOY 
corrected overlap integral are given in Table IVT1 The im- 
pact of these discussed dependencies on the 7 Bc(p,7) 8 B 
S-factor is discussed in Sect. Ivl 



IV. MOMENTUM DISTRIBUTIONS 

The stripping reaction ( 7 Be + p) + A — > 7 Bc + X 
cross section, for a specified final state of the core ( 7 Be), 
is given by [3j| (for more details, see Ref. 01) 



d 3 k r 



(2tt) 3 2Z b 



d 2 b r 



\S P (b p 



e- ik - r 5c(6 c )< A ^,^_ laIl (r)^ m6 (f) 



(8) 



where r = (r, 6*, </>) = (p,z,(j)) = r p — r c . S c (S p ) are 
the S-matrices for the core( 7 Be)+target and the pro- 
ton+target scattering, respectively. The overlap inte- 
gral qf. X is . . , , (r) describing the bound state of the 

( 7 Be + p) subsystem is defined as in Eq. (J2J) with its 
asymptotic tail corrected as described in Sect. IIIIBl 

The cross sections for the longitudinal momentum dis- 
tributions are obtained by integrating Eq. (jSJ) over the 
transverse component of k c , i.e. 



da. 



str 



dk z 



1 1 

2tt 2L -+ 



d 2 K 



l-\S p (b p )\' 



x J d 2 p\S c (b c )\ 2 

/oo 
dz exp HM <g. 6 . A _ lQ/i (r)Y hmb (f ) 
-oo 



(9) 



where k z represents the longitudinal component of k c . 

For the transverse momentum distribution in cylindri- 
cal coordinates k± = = 

dcr str 1 1 



1.2 i h.2 



d 2 k 



d 2 K 



l-\S p {b p )[ 



(2^) 2 2/ b + l J 
I dz ( d 2 p cxp(-tkj.p)S c (b c ) 



mi, 
AXJ, 



(10) 
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It is also convenient to describe the transverse momen- 
tum distributions in terms of the projection onto one of 
the Cartesian components of the transverse momentum, 
i.e. 



da. 



sir 



dk 2 



da s 



d 2 k 



(11) 



The total stripping cross section can be obtained by 
integrating either Eq. © or Eq. (|10|l . The integrals 
in Eqs. 10, (|10|l . and for the total cross sections, arc 
done numerically with use of Gaussian expansions of the 
S-matrices, as explained in Ref. |4Cj . 

The S-matrices have been obtained using the optical 
limit of the Glauber theory, i.e. 



Si (b) = exp 



kNN 



dq qpi (q) p t (q) f N N (<?) Jo (qb) 



(12) 

where pi, t (?) is the Fourier transform of the nuclear 
densities of the incident particle (i = p, 7 Be) and tar- 
get (t), and /at n (q) is the high-energy nucleon-nucleon 
scattering amplitude at forward angles, which can be 
parametrized by 

Inn (q) = ~^-vnn (i + a NN ) exp (-(3 N Nq 2 ) ■ (13) 

In this equation aNN, (Xnn, and (3nn are parameters 
which fit the high-energy nucleon-nucleon scatterin g at 
forward angles |41|. We use the values from Ref. |42| 
which have been extended to collisions at lower energies 
and corrected for isospin-average. The quantities pi (q) 
(i = p, Be) and pt (q) are calculated from radial density 
distributions taken to be of Gaussian shapes, adjusted to 
reproduce the rms radius of the proton, 7 Be, 9 Be and 
12 C, respectively. 

The Coulomb part of the eikonal phase is included ac- 
cording the prescription described in details in Ref. ^(| • 

As seen from Eqs. (0 and (|10|l . the longitudinal and 
transverse momentum distributions, as described above, 
is a direct test of the 8 B ground state wave function 
calculated with the NCSM. We do this by comparing 
to two sets of experimental data. First we compare to 
an experiment performed at MSU for the transverse mo- 
mentum distributions of 7 Be frag ments from the reaction 
8 B+ 9 Bc at 41 MeV/nucleon |4j|. In spite of the absence 
of gamma-ray coincidence data, this is a favorable case 
because the ground-state cross section dominates and be- 
cause the approximately 15% branch to the excited level 
also has 1=1 and has an almost identical shape. Be- 
cause the data is given in arbitrary units, we have multi- 
plied the contributions for the momentum distributions 
from each angular momentum channel by the same fac- 
tor so that their sum reproduces the maximum of the 
experimental distribution. Thus, by comparing our re- 
sults to this experiment we are testing the relative ratios 
between the 8 B spectroscopic factors obtained with the 
NCSM. Figures El and El shows our calculations with 
the CD-Bonn and INOY interactions, respectively The 



WS solution fit procedure was employed to correct the 
asymptotics of the NCSM overlap functions. We used 
the 10M7 model space and the optimal HO frequencies of 
h£l = 12 McV and Ml = 16 McV as determined from the 
ground-state-energy dependencies for the CD-Bonn and 
INOY potentials, respectively In both cases one notices 
that our calculations are in excellent agreement with the 
data, which they reproduce over two orders of magnitude 
in cross section. 

From Figs. ll4landlT51it is difficult to judge the quality 
of the CD-Bonn and the INOY interaction in reproducing 
the momentum distributions. Except from a noticeable 
visual difference seen for the magnitude of individual con- 
tributions for the 3 distinct 8 B states, the shape of the 
total distributions is essentially the same for the two in- 
teractions. 

We perform another comparison, this time for the lon- 
gitudinal momentum distributions, using Eq. (|3J). The 
data are taken from a GSI experiment |44| which mea- 
sured gamma rays in coincidence with 7 Be residues and 
separated the cross sections to the three final states in 
the reaction 8 B + 12 C at 936 MeV/nucleon. Figure IPol 
shows the data, which are in reasonably good agreement 
with calculations. The full drawn curves represent theo- 
retical calculations using the NCSM with the CD-Bonn 
2000 interaction. The sum of the three theoretical con- 
tributions is shown by a dashcd-linc. The experimental 
cross section is tr_ p = 94 ± 9 mb, whereas the theoretical 
one, with NCSM wave functions and the CD-Bonn 2000 
interaction, is 99.66 mb, in excellent agreement with the 
experimental one. The theoretical total momentum dis- 
tribution for this reaction is wider than the experimental 
data. The experimental FWHM (full width at half maxi- 
mum) is 95 ±5 McV/c if fitted by a single Gaussian. The 
theoretical one has a longer tail than the experimental 
one. It can be fitted by two Gaussians. The width of the 
narrower Gaussian is 96 MeV/c, also in good agreement 
with the experimental data. The reason for the longer 
tail in the theoretical momentum distribution is not well 
known, although a possible explanation is that the com- 
parison between experiment and theory, as described in 
Ref. |44|, involves a folding with the experimental reso- 
lution and a scaling to match the amplitude of the ex- 
perimental spectrum, which has not been done here. 

Figure El shows our results for exclusive p z distribu- 
tion, considering only the core excited contribution. The 
solid curve is the calculation using NCSM wave func- 
tion with the CD-Bonn 2000 interaction. The data, from 
Ref. 0, has a width of 109 ± 7 MeV/c and the theo- 
retical one is 112 MeV/c. The total experimental cross 
section for the excitation of this state is <7_ p = 12 ± 3 
mb, whereas the theoretical one is 16.36 mb, again in 
good agreement with the data. 

Figures El and El arc the same as the respective 
Figs.ElandEl but f° r NCSM wave functions calculated 
with the INOY interaction. The theoretical cross section 
for the inclusive reaction is 82.93 mb, and 13.17 mb for 
the inclusive cross section for the excited state. These 
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cross sections are smaller than the ones obtained with 
the CD-Bonn 2000 interaction. The probable reason is 
that the INOY interaction yields a smaller (? ,2 ) rms -value 
than the CD-Bonn 2000 interaction. The calculations 
yield 92 McV/c and 108 MeV/c for the for the inclusive 
and the excited state momentum distribution widths, re- 
spectively. Table IVTA summarizes the calculations for the 
total cross sections. 



V. S-FACTOR FOR 7 Bc(p, 7 ) 8 B 

We now turn our focus to the S-factor for the reaction 
7 Be(p,7) 8 B. This reaction is very important to under- 
stand the structure of our sun. The high energy neutri- 
nos from the beta-decay of 8 B come from the sun's cen- 
ter and, therefore, are a direct measure of the conditions 
in its interior (temperature, pressure, chemical composi- 
tion, etc.). The wave functions calculated with NCSM 
method, described in the previous sections will be used 
for the purpose. 

The main contribution to the S-factor for the radiative 
capture reaction 7 Be(p, 7) 8 B is due to the electric dipole 
multipolarity To calculate the S-factor, we have to 
evaluate a matrix clement of the electric dipole opera- 
tor El = ei(ri — R) with a — e for a proton and 
Cj = for a neutron and the center of mass coordinate 
R = Jr Using the notation introduced in Eqs. Q 
and (|z|) , we calculate a transition from a continuum state 

\ A ®\it£X ; Wvtj.cJiJ. M t0 thc bound statc J "> ■ 

Here, l pf c j ca i 1 .j c ( r ) is the scattering wave function of the 
7 Be+p relative motion. As the internal excitation of 7 Be 
by the electric dipole operator will have a negligible over- 
lap with the 8 B ground state, we arrive at the following 
result 



is independent on the continuous state total angular mo- 
mentum J c , <pf_j_ nT r _(r) = vF cjaaIl (r), we arrive at the 
S-factor expression 13^ 



S(E) 



47r 2 eg ff ft 2 / E + E b 
(2 J b + 1 



2(2/! + 1) ^ \-\ \ 



exp [2ttt](E)} 



jc 1 jb 



where Eb is the binding energy of the 8 B bound state 
with respect to 7 Be+p. Further, r](E) = Z a ZAe 2 /Hv p j Be 
is the Sommerfeld parameter, k is the relative momen- 
tum, E = ft 2 fc/2/i p 7Be is the center-of-mass energy of the 
proton+ 7 Be, and e, 



off = ^-^-^-i e = 3e ^ 8 ig the effec _ 

tive charge for the El capture for this specific reaction. 
The continuum wave functions are normalized so that 



r ^i c j ca iA r -> °°) - 



(16) 



where Si a j c = exp [2i<5; cJ( .], with 6i c j c (cr; c ) being the nu- 
clear (Coulomb) phase shift, and H^(r) = Gi a (r) ± 
iFi a {r), with Fi a and G; c being the regular and irregu- 
lar Coulomb wave functions, respectively. The normal- 
ization (|16|) implies that the continuum wavefunctions 
satisfy the relation 



^\^;)=5{E-E')8 V 



drr (A\J b \\EV J I l^$ (Zc i ^ ^ ; M^j^ j>) on the bound state overlap function, g 



Eq. (|15fl shows that the S-factor is a more stringent 
test of the nuclear model. Whereas the momentum dis- 
tributions discussed in the previous section depends only 

(r), 



A\J b 



Z A -l-{A-\) ATT r- 

a e VY^ 



x j[ drr a (i4AJ 6 ||ryi||*^7j; c 1 ^jM¥'w.a/ 1 j e (»-) 
_ Z A ^-(A-1) 



A 



hjb 



1 Jb Jc 



Jc Jb 



JbJc 



xj b U-l) jb+ja+h+Jo 



e -h ( jc i jb 



.1 i 
2 u 2 



X / ^^^W^^.-^) 



(14) 



with \A\Jb) the ground state of 8 B with J b = 2 + and 
\A - laii) thc ground statc of 7 Bc with If = §~. The 
7 Be proton number Za-i is equal to 4. The overlap in- 
tegral g^, X 'l\ . , (r) is defined as in Eq. J2J with its 

i'b2'^ b ' LCtli ^^^^^ 

asymptotic tail corrected as described in Sect. lIIIBl As- 
suming that the relative-motion scattering wave function 



obtained from the corrected NCSM overlap integrals dis- 
cussed in subsection IIII Bl thc S-factor also depends on 
the continuum wave function, ipfj aIi (r) . Presently, we 
have not yet developed a theory to extend the NCSM 
in order to describe continuum wave functions. There- 
fore, to obtain (pf c j caIl (r) for s and d waves, we use a 
potential model with a Woods-Saxon+Coulomb+spin- 
orbit interaction. In this respect, our approach is similar 
to that by Nollett et al. [Ifj, [46j where ah initio vari- 
ational Monte Carlo wave functions were used in com- 
bination with cluster-cluster potential model scattering 
wave functions to describe d(a,7) 6 Li, 3 H(a,7) 7 Li and 
3 He(a,7) 7 Be capture reactions. Since the largest part of 
the integrand of the last term in Eq. I15l stavs outside the 
nuclear interior, one expects that the continuum wave 
functions are well described by this model. It is possible 
to use for this purpose the same WS potential that we ob- 
tained from correcting the bound-state overlap integral. 
Then we would have a different scattering state for each 
bound-state partial wave and each NCSM model space 
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and HO frequency. In order to have the same scatter- 
ing wave function in all the calculations, we chose a WS 
potential from Ref. (3^1 that was fitted to reproduce the 
p-wave 1 + resonance in 8 B. It was argued yjj that such 
a potential is also suitable for the description of s- and d- 
waves. The parameters of this potential are summarized 
in Table EH 

We note that the S-factor results are very weakly de- 
pendent on the choice of the potential that describes the 
scattering state. This is in particular true at low ener- 
gies. Using our fitted potentials for the scattering state 
instead of the scattering potential from Table IVTI changes 
(typically increases) the S-factor by less than 1.5 eV b at 
1.6 MeV with still a smaller change at lower energies and 
no change at MeV. 

At the same time, we note that the 7 Be+p scattering 
length has been measured ^?|]- A potential model with 
no spin-orbit term was developed in Ref. 01 that fits the 
experimental scattering length. It is not consistent to use 
that potential in the present work, as we employ an al- 
ternative angular momentum coupling scheme. Still, to 
get a further insight on the S-factor energy dependence 
sensitivity to the scattering wave function, we performed 
test calculations with the s = 2 potential of Ref. • Up 
to 100 keV, we observed very little change of the S-factor 
compared to the situation when the scattering potential 
from Table IVTl was used. At higher energies, the S-factor 
increased by up to 5 eV b with a difference of 4 eV b 
at 1.6 MeV. To resolve the issue of the S-factor energy 
dependence at higher energies, we need to develop a the- 
ory that extends the NCSM to describe scattering states. 
This is a subject of current and future investigations. 

In Fig. 1201 we present our first application for the as- 
trophysical S-factor of the 7 Be(p,7) 8 B reaction. We use 
bound-state wave functions calculated with the CD-Bonn 
2000 interaction in the 10M7 model space and the opti- 
mal HO frequency Ml = 12 MeV. The WS solution fit 
procedure was employed to correct the asymptotics of 
the NCSM overlap functions. In the figure, we show the 
S-factor contributions from the dominant I = 1, j = 3/2 
and j = 1/2 partial waves by the dashed lines. Clearly, 
the j = 3/2 partial wave is the most important one. The 
upper curve, the full line, is the sum of the two contri- 
butions. The experimental data is a compilation of the 
latest experiments for the S-factor. They include direct, 
as well as some indirect measurements (Coulomb disso- 
ciation). 

The slope of the curve corresponding to the total S- 
factor follows the trend of the data. From Fig. |20| it is 
also clear that the spectroscopic factors for the 1 = 1, 
j = 3/2 and j = 1/2 partial waves are well described 
within the NCSM and the CD-Bonn 2000 interaction. 
Our calculation presented in Fig. 1201 is in a very good 
agreement with the recent direct measurement data of 
Ref. 0- We note that our S-factor energy dependence 
more resembles that obtained within the three-cluster 
model of Ref. ^(| rather than that obtained within the 
potential model of Ref. . 



In order to judge the convergence of our S-factor cal- 
culation, we performed a detailed investigation of the 
modcl-space-size and the HO frequency dependencies. 
The dependence of the 7 Be(p,7) 8 B S-factor on the size of 
the basis from 6htt to lOfiil corresponding to the above 
discussed calculation using the CD-Bonn 2000 and the 
HO frequency of Ml = 12 MeV is shown in Fig. One 
observes a very little model-space dependence with some 
oscillatory behavior. The fact that the calculation con- 
verges for a basis size within this energy interval demon- 
strates that the optimal HO frequency determined from 
the ground-state energy calculation is also close to opti- 
mal for the S-factor calculation. 

To further investigate the HO frequency sensitivity, 
we show in Fig. [221 the dependence of the 7 Be(p,7) 8 B 
S-factor on the HO frequency from Ml = 11 MeV to 
Ml = 15 MeV. The CD-Bonn 2000 NN potential and the 
1QMI model space were used. Again, the WS solution 
fit procedure was employed to correct the asymptotics 
of the NCSM overlap functions. In all cases, we obtain 
basically identical energy dependence. The absolute val- 
ues of the S-factor increase with decreasing frequency. To 
determine the optimal frequency and interpolate the con- 
verged S-factor result, we examine the basis size depen- 
dence for different HO frequencies. In Fig. [20 we show 
the Ml = 15 MeV results for model spaces from QMl 
to 10MI. We observe a steady increase of the S-factor 
with the basis size enlargement. Contrary to this situa- 
tion, the calculation using the HO frequency of Ml = 11 
MeV presented in Fig. [21 shows that the S-factor does 
not increase any more with increasing iV max . Actually, 
there is a small decrease when going from the 8M7 to 
the 1QMI. In Table IVIIII we summarize our S-factors 
at 10 keV obtained for different frequencies and model 
spaces. In addition to the results obtained using the 
WS solution fit procedure, we also present the Sn and 
the ANC obtained using the alternative direct Whittaker 
matching procedure. We note that the ANC from our ab 
initio approach are smaller than those obtained within 
the microscopic three-cluster model ^|| but still larger 
than the experimental ones from the DWBA analysis of 
Ref. (i^. In general, both procedures lead to basically 
identical energy dependence with a difference of about 1 
to 2 eV b in the S-factor with the smaller values from 
the direct Whittaker function matching procedure. Tak- 
ing into account that in the case of the direct Whittaker 
function matching the Sn increases with -/V max even at 
the HO frequency of Ml = 11 MeV, unlike in the case 
of the WS solution fit procedure, results of the two ap- 
proaches do not contradict each other. Combining all 
these results, we determine that the optimal frequency is 
between Ml = 11 and 12 MeV. Results in this frequency 
region show very weak dependence on N max , with rela- 
tive difference between the two methods always in the 
range of 5 to 8%. The full range of results is covered by 
5i 7 (10 keV) = 22.1 ± 1.0 eV b. 

We have performed a similar, although a less exten- 
sive study for the INOY NN interaction. In Fig. [2SI we 
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present the S-factor and its partial wave contributions as 
obtained using the INOY NN interaction and the ground- 
state-energy-determined optimal frequency of Ml = 16 
MeV. The \QMl model space and the WS solution fit pro- 
cedure was utilized. The energy dependence is slightly 
weaker compared to the CD-Bonn case and the S-factor 
underestimates most of the experimental data. To make 
a direct comparison with the CD-Bonn calculation, in 
Fig. |2SI we show S-factors obtained from the CD-Bonn 
2000 and the INOY NN potentials using identical HO 
frequency of Ml = 14 MeV and the 10M7 model space. 
Although the spectroscopic factors are almost the same, 
see Table the INOY S-factor is significantly smaller 
than that of the CD-Bonn 2000. This is a consequence 
of different shapes of the overlap functions displayed in 
Fig. 1131 This result is correlated with our radius and 
quadrupole moment results as well as results of momen- 
tum distributions in knockout reactions. The fact that 
the INOY NN interaction yields a smaller matter radii 
than the CD-Bonn 2000 leads to a smaller S-factor. The 
S-factor convergence for the INOY NN potential is sim- 
ilar or better than for the CD-Bonn 2000. Based on the 
results presented in the figures and the Table IVlTTl we es- 
timate the INOY S-factor result at 10 keV to be 19.0 ± 1 
eV b. 



VI. CONCLUSIONS 

We studied nuclear structure of 7 Be, 8 B and 7 ' 8 Li 
within the ab initio NCSM. We calculated overlap inte- 
grals of the 8 B ground state with 7 Be+p as a function of 
the separation between the proton and 7 Be using NCSM 
wave functions obtained in model spaces up to 10fi.fi. 
Assuming that the NCSM overlap integrals are realistic 
in the interior part, we utilized Woods-Saxon potential 
as a tool to correct their asymptotics. We performed a 
least-square fit of WS potential solutions to the NCSM 
overlap integrals in the range from fm up to about 4 
fm under the constraint that the experimental separation 
energy is reproduced. In addition, we employed an alter- 
native procedure of a direct Whittakcr function matching 



to correct the asymptotic tail of the NCSM overlap in- 
tegrals. The corrected overlap integrals were then used 
for the 7 Be(p,7) 8 B S-factor calculation. We investigated 
the dependence of the S-factor on the size of the model 
space and on the HO frequency as well as on the pro- 
cedure of the asymptotic tail correction. Based on this 
investigation, we arrived at the S-factor result at 10 keV 
to be 22.1 ± 1.0 eV b and 19.0 ± 1.0 eV b for the CD- 
Bonn 2000 and INOY NN potentials, respectively. We 
note that the CD-Bonn 2000 results for the point-proton 
rms radii, quadrupole moments and momentum distri- 
butions in knockout reactions with 8 B projectiles are in 
a better agreement with experiment than those obtained 
using the INOY NN potential. At the same time, the CD- 
Bonn 2000 predicts the point-proton rms radii of 4 He and 
6 He in agreement with experiment, while the INOY NN 
potential underpredicts both |4^. Therefore, we consider 
the CD-Bonn 2000 S-factor result more realistic. The en- 
ergy dependence as well as the zero energy value of the 
S-factor calculated using the CD-Bonn 2000 NN poten- 
tial are in a very good agreement with recent direct mea- 
surement data of Ref. Q. We stress that no adjustable 
parameters were used in our ab initio calculations of the 
8 B and 7 Be bound states. Taking into account that the 
S-factor is only weakly dependent on the potential model 
used to obtain the scattering state, we consider our re- 
sults as the first ab initio prediction of the 7 Be(p,7) 8 B 
S-factor, in particular of its normalization. 
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FIG. 1: HO frequency dependence of the 7 Be ground-state 
energy for model spaces from OhQ to 12hQ obtained using 
the CD-Bonn 2000 NN potential. The inset demonstrates 
how the values at the minima of each curve converge with 
increasing N m ax ■ 
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FIG. 2: Same as in Fig[T] but using the INOY NN potential. 
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FIG. 3: Basis size dependence of the 7 Be excitation spectrum 
in the range from OhQ to 10hQ obtained using the CD-Bonn 
2000 NN potential and the HO frequency of hQ = 12 MeV. 
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FIG. 4: Same as in Fig. 03 but using the INOY NN potential 
and the HO frequency of hil = 16 MeV. 
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FIG. 5: HO frequency dependence of the B ground-state 
energy for model spaces from Ohfl to 12fi,f2 obtained using 
the CD-Bonn 2000 NN potential. The inset demonstrates 
how the values at the minima of each curve converge with 
increasing -/V max . 




FIG. 6: Same as in Fig.QO but using the INOY NN potential. 
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FIG. 7: Basis size dependence of the B excitation spectrum 
in the range from OhSl to 10ftf2 obtained using the CD-Bonn 
2000 NN potential and the HO frequency of hQ, = 12 MeV. 
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FIG. 8: Same as in Fig. □ but using the INOY NN potential 
and the HO frequency of hil = 16 MeV. 
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FIG. 9: Overlap integral, rg(r), for the ground state of 8 B 
with the ground state of 7 Be plus proton as a function of 
separation between the 7 Be and the proton. The p-wave 
channel with j — 3/2 is shown. The full line represents the 
NCSM result obtained using the CD-Bonn NN potential in 
the lOftfi model space and the HO frequency of hfl = 12 
MeV. The dashed line represents a corrected overlap obtained 
from a Woods-Saxon potential whose parameters were fit to 
the NCSM overlap up to 4.0 fm under the constraint to re- 
produce the experimental separation energy. The dotted line 
represents a corrected overlap obtained by a direct matching 
of the NCSM overlap and the Whittaker function. 
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FIG. 10: Same as in Fig. El but for a comparison of results 
obtained using model spaces of 6, 8, 10MX Only corrections 
by the WS potential fit are shown. The fitting ranges up to 
3.6, 3.8 and 4.0 fm were used, respectively. 
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FIG. 12: Same as in Fig. M but for a comparison of results 
obtained using the HO frequencies of hQ. = 11, 13, 15 MeV. 
Only corrections by the WS potential fit are shown. The 
fitting ranges up to 4.0, 3.8 and 3.7 fm were used, respectively. 
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FIG. 13: Same as in Fig. [5] but using the HO frequency of 
hQ = 14 MeV with results obtained by the CD-Bonn 2000 
and INOY NN potentials shown for a comparison. Only cor- 
rections by the WS potential fit are shown. The fitting ranges 
up to 3.75 and 3.7 fm were used, respectively. 



FIG. 11: Same as in Fig. El but for the p-wave channel with 
3 = 1/2. 
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FIG. 15: Same as in Fig. E3] but using the INOY interaction 
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FIG. 14: Inclusive transverse-momentum distribution for 
the residue in the 9 Be( s B, 7 Be)X reaction measured at 41 
MeV/nucleon 14*3 . The theoretical calculations are based 
on Eqs. Ell and [TT] and using the CD-Bonn interaction 
for the 8 B NCSM wave function. The dotted curve is the 
sum of the individual contributions (full drawn). In or- 
der of increasing magnitude, the full curves correspond to 
I = l,j = 1/2; 7 Be(3/2-), I = l,j = 3/2; 7 Be*(l/2-), and 
I — l,j — 3/2; 7 Be(3/2 ) states, respectively. The angular 
resolution in the experiment broadens the data by approx- 
imately 4%. This has not been included in the theoretical 



data (MSU) 

■ 1 = 1, j = 1/2, 7Be(3/2-) 

■1 = 1, j = 3/2, 7Be*(1/2-) 

■1 = 1, j = 3/2, 7Be(3/2-) 
total 



INOY 




-300 -200 -100 

p y [MeV/c] 



100 



for the a B NCSM wave function. 



CD Bonn 



data (GSI) 

1 = 1, j = 1/2, 7Be(3/2-; 

= 1, j = 3/2, 7Be*(1/2-; 

= 1 , j = 3/2, 7Be(3/2-; 
total 




p z [MeV/c] 



FIG. 16: Inclusive parallel momentum distribution of 7 Be 
fragments from the reaction 8 B+ 12 C — > 7 Be+X at 936 
MeV/nucleon. The full drawn curves represent theoretical 
calculations using the NCSM with the CD-Bonn 2000 interac- 
tion. The sum of the three theoretical contributions is shown 
by a dashed-line. 
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FIG. 17: Results for exclusive p z distribution, considering 
only the core excited contribution. The solid curve is the 
calculation using NCSM wave function with CD-Bonn 2000 
interaction. 
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data (GSI) 

- 1 = 1, j = 1/2, 7Be(3/2-) 
-1 = 1, J = 3/2, 7Be*(1/2-) 
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FIG. 20: The 7 Be(p,7) 8 B S-factor obtained using the NCSM 
cluster form factors with corrected asymptotics by the WS 
solution fit. The dashed and dashed-dotted lines show the 
contribution due to the I = 1, j = 3/2 and j = 1/2 partial 
waves, respectively. The CD-Bonn 2000 NN potential, the 
lOftfi model space and the HO frequency of hi} = 12 MeV 
were used. Experimental values are from Refs. 00], El- 



FIG. 18: Same as in Fig. ED but using NCSM wave functions 
calculated with the INOY interaction. 
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FIG. 19: Same as in Fig. El but using NCSM excited state 
wave function calculated with the INOY interaction. 
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FIG. 21: The 7 Be(p,7) 8 B S-factor obtained using the NCSM 
cluster form factors with corrected asymptotics by the WS 
solution fit. The dependence on the size of the basis from 
6hQ to WHO. is shown. The CD-Bonn 2000 NN potential and 
the HO frequency of Ml = 12 MeV were used. Experimental 
values are from Refs. 0,0,13- 




NCSMhS2=ll MeV 
NCSM h£2=12 MeV 
NCSMh£i=13MeV 
NCSM h£2=14 MeV 
NCSM h£i=15 MeV 
Exp GSI-2 
ExpGSI-1 
Exp Weizmann 
Exp MSU 
Exp GSI-2 corr 
Exp Seattle 



E [MeV] 



16 



FIG. 22: The 7 Be(p,7) 8 B S-factor obtained using the NCSM 
cluster form factors with corrected asymptotics by the WS 
solution fit. The dependence on the size of the HO frequency 
from HQ. = 11 MeV to HQ = 15 MeV is shown. The CD- 
Bonn 2000 NN potential and the 10hQ model space were used. 
Experimental values are from Refs. 0,0, l|- 
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Same as in Fig. 1211 but for the HO frequency of 
MeV. 



FIG. 25: Same as in Fig. HJJbut for the INOY NN potential 
and the HO frequency of Kl = 16 MeV. 
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FIG. 26: The 7 Be(p,7) 8 B S-factor obtained using the NCSM 
cluster form factors with corrected asymptotics by the WS so- 
lution fit. The result obtained using the CD-Bonn and INOY 
NN potentials are compared. The lOhQ. model space and HO 
frequency of hi} — 14 MeV were used. Experimental values 
are from Refs. B00. 
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FIG. 24: Same as in Fig. EB but for the HO frequency of 
HQ = 11 MeV. 
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Tables 



TABLE I: The 7 Be and 7 Li ground and excited state energies 
(in MeV) obtained using the CD-Bonn 2000 and INOY NN 
potentials. The HO frequency of HQ — 12(16) MeV was used 
in the CD-Bonn 2000 (INOY) NN potential calculation. The 
ground (excited) state energies were obtained in the 12ftfi 
(lOftfi) model space. Experimental values are from Refs. l3l| . 



7 Be 





Expt. 


CD-Bonn 2000 


INOY 


i^wri)! 


37.6004(5) 


33.881 


37.161 


i i) 


0.0 


0.0 


0.0 


s *(§i f) 


0.429 


0.278 


0.501 


E *(h f) 


4.57(5) 


5.494 


5.278 


s *(f i f) 


6.73(10) 


6.999 


7.660 


S *(f 2 l) 


7.21(6) 


8.247 


8.648 


£ x (! 2 |) 


9.27(10) 


10.687 


11.331 


£x(f 2 |) 


9.9 


9.493 


10.887 


£ x (§ 2 1) 




10.120 


11.583 


s *(fi 1) 


11.01(3) 


11.717 


12.607 






''Li 






Expt. 


CD-Bonn 2000 


INOY 


|£gs(f |)l 


39.245 


35.524 


38.892 


S *(f 1 \) 


0.0 


0.0 


0.0 




0.478 


0.285 


0.513 




4.65 


5.585 


5.353 


-Mil \) 


6.60 


7.079 


7.741 



£x(§ 2 §) 


7.45 


8.522 


8.902 


S *(f 2 1) 


8.75 


9.849 


11.267 




9.09 


10.458 


11.931 


1) 


9.57 


11.033 


11.691 


S *(fl I) 


11.24 


11.974 


12.826 



TABLE II: The 7 Be and 7 Li point-proton rms radii (in fm), 
ground-state quadrupole (in efm 2 ) and magnetic (in /xn) mo- 
ments and Ml transitions (in fj^) obtained within the NCSM 
for different HO frequencies (given in MeV) and model spaces 
for the CD-Bonn 2000 and INOY NN potentials. Experimen- 
tal values are from Refs. pjlll . 

7 Be 

CD-Bonn 2000 



HQ iV ma x 


r P 


Q 




B(M1; i~ -»• 


I") 


12 6 


2.311 


-4.755 


-1.150 


3.192 




12 8 


2.324 


-4.975 


-1.151 


3.145 




in in 


z.o4z 


-5.153 


-1.141 


9 1 1 A 
0.114 




12 12 


2.365 










11 6 


2.377 


-5.029 


-1.150 


3.203 




11 8 


2.377 


-5.186 


-1.151 


3.162 




11 10 


2.383 


-5.354 


-1.155 


3.125 




Expt. 


2.36(2) 


-1.398(15) 
INOY 


3.71(48) 




HQ iV ma x 


r P 


Q 




B(M1; I" 


?, > 


16 6 


2.114 


-3.946 


-1.161 


3.158 




16 8 


2.149 


-4.212 


-1.157 


3.119 




16 10 


2.181 


-4.459 


-1.151 


3.092 




16 12 


2.214 










14 6 


2.190 


-4.245 


-1.162 


3.180 




14 8 


2.209 


-4.454 


-1.158 


3.139 




14 10 


2.227 


-4.654 


-1.155 


3.108 




Expt. 


2.36(2) 




-1.398(15) 


3.71(48) 
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CD-Bonn 2000 






HQ, iV ma x 


r P 


Q 




B(M1; \" 


2 ") 
2 > 


12 6 


2.149 


-2.717 


+3.027 


4.256 




12 8 


2.156 


-2.866 


+3.020 


4.188 




12 10 


2.168 


-3.001 


+3.011 


4.132 




12 12 


2.188 


-3.130 








Expt. 


2.27(2) 


-4.06(8) +3.256 
INOY 


4.92(25) 




hQ, iV max 


r P 


Q 




B(M1; I" 


?, > 


16 6 


1.963 


-2.334 


+3.039 


4.192 




16 8 


1.990 


-2.500 


+3.029 


4.136 




16 10 


2.015 


-2.648 


+3.021 


4.098 




16 12 


2.042 


-2.788 








Expt. 


2.27(2) 


-4.06(8) 


+3.256 


4.92(25) 





TABLE III: The B and Li ground and excited state energies 
(in MeV) obtained using the CD-Bonn 2000 and INOY NN 
potentials. The HO frequency of hQ. — 12(16) MeV was used 
in the CD-Bonn 2000 (INOY) NN potential calculation. The 
ground (excited) state energies were obtained in the 12HQ 
(10ftf2) model space. Experimental values are from Refs. [32^ . 



8 B 

Expt. CD-Bonn 2000 INOY 



17 



|£ gs (2+l)| 


37.7378(11) 


32.284 


36.148 


£ x (2+l) 


0.0 


0.0 


0.0 


E K (lfl) 


0.774(6) 


0.804 


1.199 


£ x (3 + 1) 


2.32(20) 


2.977 


2.854 


£ x (0+l) 




2.229 


3.853 


Sx(l|l) 




2.988 


4.540 


S x (2+1) 




3.824 


4.897 


S x (l+1) 




4.827 


6.459 


S x (2+1) 




5.175 


5.908 


-Ex (4+1) 




6.482 


7.138 


S x (3+1) 




7.325 


8.572 


£x(0+2) 


10.619(9) 


10.782 


11.926 




8 Li 








Expt. 


CD-Bonn 2000 


INOY 


|£ gs (2+l)| 


41.277 


35.820 


39.938 


S x (2+1) 


0.0 


0.0 


0.0 


^x(l + l) 


0.981 


0.855 


1.264 


E x (3+1) 


2.255(3) 


3.019 


2.871 


S x (0+1) 




2.480 


4.225 


Sx(l+1) 


3.21 


3.247 


4.903 


S x (2+1) 




3.977 


5.114 


£ x (l+l) 




5.023 


6.758 


£ x (2+l) 




5.290 


6.071 


E K (4+1) 


6.53(20) 


6.691 


7.398 


£ x (3+l) 




7.570 


8.915 


S x (0+2) 


10.822 


10.898 


12.049 



TABLE IV: The 8 B and 8 Li point-proton rms radii (in fm), 
ground-state quadrupole (in efm 2 ) and magnetic (in /in) mo- 
ments and Ml transitions (in fj,^) obtained within the NCSM 
for different HO frequencies (given in MeV) and model spaces 
for the CD-Bonn 2000 and INOY NN potentials. Experimen- 
tal values are from Refs. 1321 . 
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B(M1; 1+ ->■ 
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2.436 


+5.218 


+1.463 


3.498 




12 
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2.463 


+5.420 


+1.455 


3.506 




12 


10 


2.487 


+5.636 


+1.455 


3.490 




12 


12 


2.520 










11 


6 


2.514 


+5.525 


+1.515 


3.491 




11 


8 


2.528 


+5.696 


+1.501 


3.495 




11 


10 


2.542 


+5.871 


+1.496 


3.539 
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Q 




B(M1; 1+ ->■ 
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16 
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2.199 


+4.049 


+1.192 


3.669 




16 


8 


2.241 


+4.306 


+1.207 


3.669 




16 


10 


2.277 


+4.580 


+ 1.227 


3.649 




16 


12 


2.317 










15 


6 


2.241 


+4.242 


+1.238 


3.681 




15 


8 


2.276 


+4.468 


+ 1.244 


3.684 




15 


10 


2.305 


+4.710 


+ 1.257 


3.667 




14 


6 


2.291 


+4.467 


+1.288 


3.696 




14 


8 


2.318 


+4.660 


+1.286 


3.702 




14 


10 


2.340 


+4.867 


+1.291 


3.688 
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B(M1; 1+ -> 


2+) 



18 



12 6 


2.139 


+2.588 


+1.238 4.454 


12 8 


2.139 


+2.690 


+1.243 4.428 


19 in 

1Z 1U 


9 1/1^ 

Z. l^iO 


I O 'TO A 
-T~Z. / 


9/11 A *}Q*} 
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2.161 






Expt. 


2.26(2) 


+3.27(6) 


+1.654 5.01(1.61) 






INOY 




HQ, iVnmx 


r P 


Q 


y, B(M1; 1+ -> 2+) 


16 6 


1.938 


+2.279 


+1.469 4.635 


16 8 


1.956 


+2.377 


+1.456 4.610 


16 10 


1.972 


+2.477 


+1.439 4.578 


16 12 


1.991 






Expt. 


2.26(2) 


+3.27(6) 


+1.654 5.01(1.61) 


TABLE V: 


The < 8 B(2+)| 7 Be(/D 


+p(l,j)) spectroscopic fac- 


tors obtained within the NCSM for different HO frequencies 


and model 


spaces for the CD-Bonn 2000 and INOY NN po- 


tentials. 








CD-Bonn 2000 




7 Be+p 


hQ, [MeV] 


^ max 


F(i,§) 


9 \ A ! 9 I 2 V- 1 -' 2 / 


11 


6 


0.977 


0.120 0.285 


11 


8 


0.967 


0.116 0.280 


11 


10 


0.959 


0.111 0.275 


12 


6 


0.978 


0.107 0.287 


12 


8 


0.969 


0.103 0.281 


12 


10 


0.960 


0.102 0.276 


14 


6 


0.979 


0.086 0.288 


14 


8 


0.967 


0.085 0.284 


14 


10 


0.958 


0.085 0.280 


INOY 




'Be+p If 


hQ, [MeV] 


^ max 


3-Q 1) 

9 \ ' 11 


3"(1 I) l-fl 3n 

9 \ ' 2 J 9 \ ' 9' 


14 


6 


0.987 


0.074 0.283 


14 


8 


0.976 


0.072 0.277 


14 


10 


0.966 


0.072 0.272 


16 


6 


0.988 


0.060 0.281 


16 


8 


0.977 


0.061 0.276 


16 


10 


0.965 


0.063 0.271 



TABLE VI: Parameters of the Woods-Saxon potentials 
obtained in the fits to the interior part of the NCSM 
( 8 B(2 + 3 )| 7 Be(/ 1 r )+p(Z, j))(r) overlap functions under the con- 
straint to reproduce experimental separation energies. The 
p-wave channels for the 7 Be ground and the first excited state 
are shown. Results for the Whfl model space and HO frequen- 
cies of HQ, = 12(14) MeV for the CD-Bonn 2000 (INOY) NN 
potential. Woods-Saxon potential parameters from Ref. l3Sl 
(with slight modification of Vq) that we use for the scattering 
states are also shown. 

CD-Bonn 2000 lOHQ hQ. = 12 MeV 

I? (I J) Vo Ro "Q Yh RlB aiB Rc 

|"(1,|) -51.037 2.198 0.602 -9.719 2.964 0.279 2.198 

§~(1,±) -45.406 2.613 0.631 -8.414 2.243 0.366 2.613 

|~(1, |) -49.814 2.235 0.553 -17.024 3.080 0.338 2.235 
INOY WHO, HQ. = 14 MeV 

Vq Rq a V la R la a la Rg 

|"(1,|) -58.836 2.052 0.561 -7.518 2.768 0.253 2.052 



3 
2 


"(1, 


\) 


-55.924 


2.470 0.580 


-17.454 


2.027 


0.429 


2.470 


1 
2 


"(1, 


1) 


-44.300 


2.455 0.509 


-13.325 


1.011 


0.347 


2.455 










scattering 


; state 














V 


Ro ao 


Ms 


Ris 


ais 


Rc 








-42.2 


2.391 0.52 


-9.244 


2.391 


0.52 


2.391 



TABLE VII: Cross sections for the proton-removal reactions 
8 B + 9 Be at 41 MeV/nucleon (MSU) and 8 B + 12 C at 936 
MeV/nucleon (GSI). The calculated total inclusive cross sec- 
tions are given by a^. for the CD-Bonn 2000 interaction 
and Cjnc for the INOY interaction. The lOftfi model space 
and HO frequencies of 12(16) MeV were employed for the 
CD-Bonn 2000 (INOY). The calculated cross sections for the 



excited state are given by and <7exc , respectively. 





(exp) 
inc 

[mb] 


(th-B) 
Vine 

[mb] 


(th-i) 

inc 

[mb] 


(exp) 
O exc 

[mb] 


(th-B) 

" exc 

[mb] 


(th-i) 

t-* exc 

[mb] 


MSU 




82.96 


71.85 




15.31 


13.26 


GSI 


94 + 9 


99.66 


82.93 


12 + 3 


16.36 


13.17 



TABLE VIII: The calculated 7 Be(p,7) 8 B S-factor, in eV b, at 
the energy of 10 keV. Two ways of correcting the NCSM over- 
lap asymptotics, by the Woods-Saxon potential solution fit 
(WS) and by a direct Whittaker function matching (Whitt), 
are compared. The asymptotic normalization constants, in 
fm _1//2 , correspond to the Whittaker function matching case. 
Results obtained using the CD-Bonn 2000 and INOY NN po- 
tentials at different HO frequencies and model spaces as well 
as the NCSM extrapolated Svt values with their estimated 
errors are presented. 



CD-Bonn 2000 



hQ, [MeV] 


N max 


Cl,3/2 


Cl.l/2 


cjWhitt 
<->17 


qWS 


15 


6 


0.647 


0.195 


16.81 


17.80 


15 


8 


0.660 


0.206 


17.58 


18.87 


15 


10 


0.672 


0.216 


18.33 


19.81 


14 


10 


0.680 


0.220 


18.78 


20.21 


13 


10 


0.692 


0.234 


19.64 


21.02 


12 


6 


0.693 


0.240 


19.75 


21.24 


12 


8 


0.696 


0.242 


19.96 


21.14 


12 


10 


0.704 


0.247 


20.45 


21.66 


11 


6 


0.715 


0.261 


21.30 


22.38 


11 


8 


0.715 


0.263 


21.33 


23.04 


11 


10 


0.720 


0.262 


21.60 


23.06 




NCSM 517(10 keV) 




22.1 ± 1.0 






INOY 








Tifi [MeV] 


-^V max 


Cl,3/2 


Cl.l/2 


a Whitt 
£>17 


qWS 


16 


10 


0.641 


0.182 


16.34 


17.49 


15 


10 


0.649 


0.189 


16.83 


17.95 


14 


6 


0.652 


0.190 


16.94 


17.78 


14 


8 


0.654 


0.194 


17.12 


18.21 


14 


10 


0.660 


0.198 


17.44 


18.46 




NCSM 5*17(10 keV) 




19.0 + 1.0 



